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ABSTRACT 

Measurements of temporal and spatial scale sizes of small scale 
temperature fluctuations were made in 19 meters of water from the 
NUC Oceanographic Research Tower off San Diego, California. Hor- 
izontal and vertical measurements were made using a two meter 

straight line array of seven thermistors. In addition to temperature, 

/ 

waves and orthogonal water particle velocities (turbulence) were meas- 
ured simultaneously. Internal waves were found on two of five runs. 

Average correlation length and time computed for the horizontal 
array runs without internal waves were 118 cm (standard deviation 
21.5 cm) and 10 seconds (standard deviation 4 seconds). Measure- 
ments suggested that correlations obtained from the mixed layer were 
longer than those obtained at the edge of the thermocline. Internal 
wave induced turbulence reduced spatial and temporal scales to 33 cm 
and 4 seconds (standard deviation 0.7 seconds). 

Vertical measurements taken in the presence of internal waves 
yielded spatial and temporal correlations of 38 cm and 4 seconds 
(standard deviation 0.7 seconds). 
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I. INTRODUCTION 



A. HISTORICAL 

For the past two years research has been carried out at the Naval 
Postgraduate School Department of Oceanography to study the shallow 
water interactions of various oceanographic processes that occur at 
small scales. Prior to the conduct of this experiment in May 1973, 
experiments were made in October 1971 and June 1972 from the Naval 
Undersea Center (NUC) Oceanographic Research Tower off Mission 
Beach, San Diego, California. Previous experiments included the 
simultaneous measurements of several parameters such as surface 
waves, water particle velocities, and fluctuations of temperature, 
sound speed, sound phase, and salinity. The results of these exper- 
iments are presented for the most part in MS theses of Naval Post- 
graduate School students. Smith [1971], Rautman [1971], Fitzgerald 
[1972], and Alexander [1972] reported on sound amplitude and phase 
fluctuations. Bordy [1972] made wave and water particle motion 
measurements. Krapohl [1972] made more detailed wave and turb- 
ulence measurements. Seymour [1972] initiated work on salinity 
measurements which were completed by Frigge [1973]. Duchock [1972] 
attempted to correlate sound velocity measurements with temperature 
and other fluctuations. Haley [1973] re-examined the October 1971 
experiment in an attempt to interrelate all parameters. Gossner [1973] 
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compared measured sound speed at several depths with sound speed 
computed from empirical equations. Whittemore [1973] made pro- 
gress in relating the near surface temperature field to wave motion. 

From these studies, it became apparent that a more detailed 
examination of the small scale temperature field was required before 
further progress could be made in determining the interactions of 
other parameters. 

B. STATEMENT OF THE PROBLEM 

The creation of the small scale temperature field in shallow water 
is a complex process and may occur at various depths in the water 
column. At the surface, wind stresses create waves and if the stress 
is sufficient to produce wave breaking (white caps), then a turbulent 
transfer of momentum downward will result with an attendant down- 
ward flux of any excess heat in the surface water. Wind stress will 
also produce surface drifts. The shears from these drifts produce a 
certain amount of turbulence which may also be transmitted downward. 
This process of surface originated turbulent mixing extends downward 
to a depth proportional to the wavelength of the surface waves. 

The ocean bottom also exerts a stress on water moving over it 
due to, for instance, tidal currents. The resulting shear zone causes 
turbulent mixing which affects the small scale temperature field and 
extends upward a distance proportional to the stress exerted. 
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The presence of a thermocline inhibits the vertical extent of the 
turbulence caused by the processes described above. Since the degree 
of stability of the water column tends to define the degree of vertical 
motion, the increase in stability implied by the density gradient will 
reduce the vertical motion of the turbulence and thus the vertical 
mixing. 

It is paradoxical that the thermocline, which inhibits the turbulent 
mixing due to surface and bottom shears, may itself be the point of 
creation of turbulence. Internal waves propagate along the thermocline 
and can lead to the creation of intense turbulence. If the water masses 
separated by the thermocline experience relative motion, then the shear 
across the thermocline will produce turbulence. This relative motion 
may be caused by tidal currents or wind drifts. [La Fond, 1962] 

It is interesting that turbulent mixing from surface and bottom 
shears tends to intensify the thermocline whereas turbulence caused 
by internal waves or current shears tends to dissipate it. Wind 
induced mixing redistributes heat downward to the thermocline so 
the deep water becomes warmer; bottom mixing essentially re- 
distributes heat downward away from the thermocline so the upper 
level becomes cooler. The gradient thus becomes stronger. The 
internal waves tend to weaken the gradient through complex non-linear 
wave interactions and make the temperature difference smaller. 

The greater the shear across the thermocline, the greater the 
possibility of instability and turbulence. It is well known that the 
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intensity of the shear determines the rate of growth of an unstable 
disturbance in a steady, stable, parallel, shearing motion. The shear 
instability may be described by the local Richardson number, J(z), 




3/f. Hit 

C^) L 



where g is the acceleration due to gravity, £ is the reference 
dens ity, Hit is the rate of density change across the thermo- 

cline, and Viz is the velocity shear across the thermocline. When 
the Richardson number is less than a critical value (J(z)^ 0.25) then 
turbulence is possible though not certain [Phillips, 1969]. Since the 
maximum shear is normally found at the crests and troughs of waves, 
this turbulence should be expected to occur most of the time or with 
greatest intensity at these points. If the thermocline is intense, the 
critical shear will be reached earlier than if it were diffuse; an 
internal wave in an intense thermocline will dissipate its energy into 
turbulence more quickly than one in a diffuse thermocline. The tur- 
bulent mixing resulting from the passage of an internal wave results 
in the reduction of the shear in the thermocline to a value less than 
the critical value. Creation of turbulence by this means will then 
presumably cease until the thermocline is again restored. 

An internal wave need not be a turbulent phenomenon. If there is 
no turbulence being created at the thermocline, then the regenerating 
forces of turbulent mixing from surface and bottom shears may 
restore the intensity again. 
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The interplay between propagation of energy in an internal wave 
and its dissipation may be summarized in the following way. An 
internal wave is initiated by some process and propagates along a 
density gradient. The shear at the boundary increases possibly due 
to the wave entering shallow water and changing its form or possibly 
due to an increase in the intensity of the density gradient or thermo- 
cline. If this shear exceeds a critical value, then turbulence may be 
created, probably at the crest or trough. The resulting turbulent 
mixing reduces the gradient possibly below the critical value. The 
wave then either continues propagating with energy reduced by losses 
to turbulence, or, if the shear is not reduced below the critical value, 
the wave is entirely dissipated. Areas of density gradients reduced 
by the mixing of past waves will not be able to support waves of such 
intensity of shear. 

Theoretically, then, there is a cycle between an intense thermo- 
cline with intense internal wave activity producing turbulence, and a 
diffuse thermocline with internal waves producing little or no 
turbulence. 

The kinetic energy of the turbulence is eventually dissipated by 
entrainment. Shears at the boundaries of the turbulent eddies formed 
by the mixing processes described above cause the shedding of smaller 
eddies, among which the energy of the original eddy is distributed. 

The energy is cascaded through eddies of decreasing size and energy 
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until eventually the energy is no longer sufficient to overcome the 
viscous forces of the surrounding water. The remaining energy is 
then dissipated as heat energy. 

Small scale temperature gradients may exist in a field even 
though all velocity gradients have been removed by viscosity because 
the scales at which molecular diffusion of heat becomes effective are 
much smaller than those at which molecular viscosity is effective. 
This small scale temperature field is subject to being moved about 
as a passive scalar riding on the water particle motion caused by 
surface and internal waves. 

C. OBJECTIVES 

An experiment was conducted in May 1973 in order to make a 
careful examination of the small scale temperature field in the near 
surface region where the influence of surface waves is felt. The 
principal sensors were a moveable array of thermistors, a wave 
gauge, and a flowmeter. The object was to verify Whittemore's 
suggested wave- temperature coupling mechanism. 

Particular attention was given to determining time and length 
scales. The scales that can be resolved are determined by the char- 
acteristics of the sensors and their spatial relationships with each 
other. For this experiment spatial scales were determined by the 
minimum thermistor spacing (about eight centimeters) and by the 
maximum spacing (about two meters). Minimum temporal scales 
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depended on the response time of the thermistors (0. 015 seconds). 
Maximum scales depended on the length of the record. A maximum 
scale of five percent of the record length of about 20 minutes was used 
for the maximum lag time in calculating correlation functions giving 
an upper limit of one minute period waves. 

The conditions of the near surface region were fulfilled by the 
shallow water depth of approximately 19 meters. 

The May date was chosen partly for logistic reasons but also 
because only moderate wave activity was expected and vertical tem- 
perature gradients would not be intense. Internal wave activity was 
not expected to be a major source of turbulence. Unfortunately, 
conditions did not occur as predicted. Wind waves and swell were 
small or nonexistent and internal waves dominated the temperature 
scale much of the time. Although these conditions detracted from the 
original aim of the experiment, they provided an excellent opportunity 
to examine small scale (turbulent) temperature fluctuations in the 
presence of internal waves. 

This thesis presents an analysis and discussion of the data taken 
during the May 1973 experiment conducted from the NUC tower. 
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II. DESCRIPTION OF THE EXPERIMENT 



A. LOCATION 

The experiment was conducted from the Naval Oceanographic 
Research Tower operated by the Naval Undersea Center (NUC) 
located about one mile off Mission Beach, San Diego, California. 

The tower consists of four deck levels about 40 ft square each on a 
frame of four legs. One of the decks contains a small deck house 
suitable for housing electronic equipment. Vertically mounted rails 
from one of the lower decks run down three sides of the support frame; 
equipment may be mounted on these rails and lowered to any pre- 
selected depth. The rails on the seaward face were used for the 
experiment since wave action was expected to come from that 
direction. The tower offers stable support and ample working space. 
Water depth at the tower is about 18. 3 m with a smooth sandy bottom. 
[La Fond, 1965] 

B. TIME 

The experiment was conducted 15 and 16 May 1973; this time of 
year is a time of transition between winter and summer conditions. 

Air and sea surface temperatures are normally in the vicinity of 
16. 7°C, low wind, sea, and swell activity is usual, the thermocline 
is not likely to be strong, and internal wave activity is not yet intense. 
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Actual meteorological and oceanographic environmental factors 



observed were: 



Day 1: 



RH: 

SST: 



air temp: 
dew point: 




overcast 
wind: calm 

sea sfc: negligible wind 
waves; swell 0. 7 m from west 



16.5 C 



Day 2: 



RH: 

SST: 



air temp: 
dew point: 



14.4-15.5 C 
13. 3-13. 9°C 
84% 




overcast 
wind: calm 

sea sfc: capillary waves; 

swell 0.6 m from west 



Conditions were essentially constant over the two days except the sky 
cleared at the end of the last run on the second day. Nansen casts 
were made for all runs, results are in table I; bathythermograph and 
sea surface temperature measurements were made for each run with a 
shallow BT and a bucket thermometer. The calibration of the BT was 
uncertain and there was considerable hysteresis. However, an 
estimation of thermal gradients was made for each run; the results 
are in table II and figure 1. 

Tidal conditions were as follows: 



Day 1 



1340 low 
1955 high 



1125 ebb 
1729 flood 



1407 slack 



Day 2 



0917 high 
1409 low 



0933 slack 
1432 slack 



1158 ebb 
1754 flood 
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TABLE I 



RESULTS OF NANSEN CASTS 



RUN 



NUMBER 


TIME 


DEPTH 


SALINITY 




1 


1315 


3 m 


33.743 0/00 




2 


1440 


8 m 


33. 736 0/00 




3 


1600 


14 m 


33.739 0/00 




4 


1730 


3 m 


33. 738 0/00 




6 


1100 


3 m 


33. 743 0/00 








TABLE II 








ESTIMATION OF THERMAL 


GRADIENTS 




RUN 




TOP OF 


INTENSITY OF 


NUMBER 




GRADIENT 


GRADIENT 


1 




4. 25 m 


o 

0. 98 


C/m 


2 




10. 5 m 


o 

0.88 


C/m 


3 




7.6 m 


o 

0. 79 


C/m 


4 




14. 7 m 


0.25° 


C/m 


6 




sfc 


0.44° 


C/m 
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Temperature (°C) 




Temperature (°C) 




( I indicates depth of frame) 



Figure 1. Bathythermograph traces 
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Depth (m) Depth (m) 

18.3 0 18.3 



C. EQUIPMENT USED 



A specially constructed 6x6x8 foot frame was mounted on the 

vertical rails of the seaward face of the tower by means of sliding 

shoes. An array of seven individually movable thermistors and a two 

channel electromagnetic flowmeter were mounted on the outside face 

of the frame. The frame was then lowered to the predetermined depth 

by means of an electric winch mounted on the deck at the level of the 

railhead. Approximate depth was determined by means of a marked 

cord attached to the top of the frame and payed out from the winch 

deck. An I-beam was mounted horizontally on the seaward face of 

the deck house extending over the equipment frame. A trolley ran on 

the underside of the beam; a Baylor wave gauge was suspended between 

the trolley and the top of the outside face of the equipment frame 

[See Fig. 2]. S ince the frame was displaced horizontally as it slid 

o 

down the rails due to a 5 cant in the rails, the trolley was moved to 
keep the wave gauge vertical and directly over the flowmeter and ther- 
mistor array. Cable runs ran from all equipments to the recording 
room in the deck house. 

The core of the experiment was seven thermistors. The ther- 
mistors were type K496 from Fenwall Electronics Inc. , designed for 
oceanographic work. The thermistors consist of two specially aged 
beads in series; since there is a statistically small probability of both 
beads drifting or drifting in the same direction with time, they are 
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I-beam 




Figure 2. Layout of the Experiment 
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exceptionally stable and have a high repeatable accuracy. The thin 
glass watertight sealant around the beads allows a short response 
factor of 0.015 sec. High resolution is obtained from the high resist- 
ance (15 K ^ at 25°C) and a high resistance to temperature change 
o o 

ratio (540 K ft / C at 25 C). Due to low self heating, measurements 
accurate to the nearest 0.01°C can be obtained in seawater [Shonting, 
1968]. The thermistors were mounted in a 0. 3x9 cm aluminum tube 
with the beads extending just beyond the lip of the tube. The tube was 
then set in a 6 cm diameter epoxy base in which was sealed the con- 
nection between the cable run and the thermistor leads. The beads 
were protected in transit by a stainless steel sleeve that fitted over the 
tube; during the experiment they were protected from fish and seaweed 
by a 3 mm diameter bronze cage (Plate 1). Pipe straps were attached 
to the epoxy bases which allowed the thermistors to be mounted on the 
frame. 

Associated with each thermistor was a Wheatstone bridge. Of the 
four legs of the bridge, two were occupied by equal fixed resistors; 
the balancing leg consisted of a 1 turn and a 10 turn potentiometer 
for coarse and fine tuning; the thermistor formed the fourth leg 
(Fig. 3). The bridge was constructed so that it could be balanced 
with the thermistor in place. A decade box could then be substituted 
for the thermistor and the bridge rebalanced by means of the decade 
box. The balancing resistance and thence the temperature of the 
thermistor would then be known. With the thermistor back in place, 
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m 



Plate 1 

Thermistor and Guard 



I ■ 



-s-;l 



m2 



m ri 



5 



cm 
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Battery# 6.5 V 




Figure 3. Circuit Diagram for Wheatstone Bridges 
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Calibrate 



any voltage fluctuation across the bridge would be a function of the 
thermistor temperature. Power was supplied by a 6. 5 V mercury 
battery. 

The thermistors were calibrated by immersion in a fresh water 

bath of known temperature. The temperature of the bath was measured 

by a quartz thermometer which had itself been calibrated at 0°C in a 

distilled water and ice bath. Plots were made of bridge output versus 

o o 

temperature for temperatures between 10 C and 21 C in one degree 
increments. The curves were sufficiently linear for a valid assump- 
tion of linearity over a range of a few degrees (Fig. 4). 

The calibration factor for the individual thermistors was then 
determined by solution of the Wheatstone bridge equations: 



E=(R D +R 1 )i 2 



E=(R T +R 2 )i 1 




or 



V= r d e 

(R d + r i 




R'pE 



where R^, is the instantaneous resistance of the thermistor, R^ is the 



balance point of the bridge as determined by the decade box, R^ and 
R 2 are the constant resistances, and E is the applied voltage, and V 



is the output of the bridge when unbalanced. The calibration factor 
thus determined was typically 14. 9°C/V, (Fig. 5). 



Resistance (Kil) 




Temperature (°C) 

Figure 4. Typical Temperature vs 
Resistance Calibration Curve for Thermistors 
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Figure 5. Calibration Factor for Wheatstone Bridge 
Balanced at 14°C 
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The Baylor wave gauge used has an accuracy of within 0. 01 of 
wave height and a response time of less than 0. 06 sec. The gauge 
used was 4. 1 m long; its calibration factor was 7. 7 m/V. The cal- 
ibration curve is figure 6. 

The two channel electromagnetic flowmeter was an Engineering 
Physics Company Water Current Meter model EMCM3B. The ac- 
curacy is to within 0.015 m/sec and has a flat response out to at 
least 0.5 seconds. The calibration factor was 0.541 m/sec/ V. A 
complete description of this flow meter and calibration is given by 
Steer [1972]. 

D. RECORDING 

All signals except those from the thermistors were recorded 
directly from the instruments. The Wheatstone bridges had to be 
balanced with the thermistors in the ocean during the procedures 
described for calibration. Voltage fluctuations caused by the tem- 
perature fluctuations in the ocean were never serious enough to 
prevent balancing the bridges. 

Since the thermistor signals were very small (on the order of 
millivolts), they had to be amplified before being recorded. Seven 
identical amplifiers were not available, so three Preston 8300XWB 
and four Hewlett-Packard 2470B amplifiers were used. Unfortunately, 
the only gain setting they had in common was xlO. For some runs in 
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Immersion Length (m) 




Figure 6. Calibration Curve for Baylor Wave Gauge 
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water with little thermal activity a higher gain setting could be used: 
x30 on the Preston amplifiers and x50 on the Hewiett-Packards. In 
the presence of internal waves when temperature excursions were 
relatively great, the lower gain setting was used to avoid overdriving 
the tape recorder. 

The Preston amplifiers were equipped with a filter; consequently, 
on the first day channels with the Preston amplifiers were filtered at 
10 Hz low pass to filter out any 60 Hz noise. However, on the second 
day, the filters were not used in order not to phase shift the signals. 

All signals were then recorded simultaneously on a Sangamo 
3562 14 channel FM magnetic tape recorder. Full scale deflection 
was set at + 1 V rms. Tape speed used was 1 7/8 ips; this allowed 
maximum speedup (to 60 ips) during playback. 

E. EXPERIMENTAL PROCEDURE 

Equipment set up for the two days was similar in most respects. 

On the first day the thermistor array was vertical as shown in figure 7. 
On the second day the array was horizontal for the first two runs 
(Fig. 8), then for the last three runs, thermistor one was moved 
below the remainder of the horizontal array to give some indication 
of vertical temperature gradient. 
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#4 
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1 
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i 

1 




#5 
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1 

4/ 


Flowmeter 


8.9 cm 

'■4 > 








A 

#6 


V 

In 

o 

3 

Thermistors 






# 7 


• 162 . 5 cm 



Figure 7. Configuration of Array for 
Runs One, Two, Three, and Four 
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Figure 8. Configuration of Array for Run Six 
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162.5 



Run times and depth of the top of the frame are given below: 



Run number 

1 

2 

3 

4 

5 

6 

7 

8 
9 



Time 

1300-1340 

1437-1510 

1607-1640 

1722-1755 

1017-1034 

1050-1201 

1244-1315 

1339-1400 

1410-1430 



Depth (m) 

2.7 
5.2 
11.3 

2.7 

2.7 

2.7 
8.5 
13 . 1 
13 . 1 
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HI. DATA ANALYSIS 



A. DATA PREPARATION 

Analysis of the data was done at the Naval Postgraduate School. 
Initially all signals were graphically displayed on a strip chart for 
visual inspection, polarity checks, and checks for errors or anomalies. 
Immediately, two runs were discarded. Runs five and eight had been 

started with amplifier gains too high and too low and had overdriven 

» 

and underdriven the tape recorder respectively. Runs six and nine 
had been made at the same depths and with the same configuration and 
therefore replaced five and eight. On runs seven and nine, there was 
a great deal of noise on the signal from the sixth thermistor; these 
runs were not analyzed. 

Some long term fluctuations were noticed on all runs. Two runs 
showed significant internal wave activity. Run three showed a distinct 
internal wave in the top three thermistor traces (Fig. 9)* Run six 
was started in the middle of an internal wave, spanned a second distinct 
wave, and included two patches of increased turbulence that may have 
been the remains of other internal waves (Figs. 10, 11). 

Twenty minute portions of each run were chosen for analysis. 

This was the first 20 minutes in all but run six. In run six the portion 
for analysis did not include the large internal waves; the microstructure 
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35 





36 




37 




in the internal waves did not appear to be representative of the rest 
of the run. 

Once suitable portions of the data had been identified, they were 
digitized on a hybrid computer CI5000, XDS9300. Some runs required 
filtering of 60 Hz noise; this was done with Krohn-hite 3750 filters. 
Playback was at 60 ips with a sampling rate of 400 Hz, corresponding 
to a real time sampling rate of 12.5 Hz. The signals were further 
amplified during digitization to take advantage of the full dynamic 
range of the computer and thus increase resolution. Amplification was 
either x50 or xlOO depending on the amount of fluctuation in the signal. 

The seven track digital tape thus made on the hybrid computer 
was then transformed onto a nine track digital tape compatible with 
the school's IBM 360-67 computer on which the analysis was 
accomplished. 

Since the scales to be analyzed were on the order of about one to 
100 seconds, it was decided to filter out digitally any activity with 
a period of 100 seconds or greater. A Fourier analysis was made 
of each 20 minute segment of data; the coefficients were used to 
generate a time series describing the waves of period 100 seconds and 
greater. This series was then subtracted from the original time series. 
It was the filtered time series that was used for analysis. 

The DC component and any trend were then deleted. 
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B. METHOD OF ANALYSIS 

Temporal auto- and cross-correlations were computed for selected 
sensors, A Parzen window was first applied to the correlation func- 
tions. Auto- and cross-spectra were computed by a Fourier trans- 
formation on the correlation functions. Coherence and phase were 
computed from the cross-spectra. Due to the desire to cut down on 
computation time but keep the 20 minute record length, the number of 
data points was arbitrarily cut in half by deleting every other data 
point. This reduced the sampling rate to 6.25 Hz resulting in a 
Nyquist frequency of 3. 1 Hz. The reduction in the Nyquist frequency 
was not considered critical, since the spectra became essentially 
flat and the signal lost in noise before the 3. 1 Hz frequency was 
reached. Aliasing was experienced possibly from 60 Hz noise which 
was never effectively filtered out. The energy was always easy to 
identify, however, and appeared in narrow spikes at 2. 59 Hz, and 
occasionally at 2. 09, 1. 04, and 0. 67 Hz. 

Spatial correlation was computed by averaging cross variance 

(C. V. ) between the first thermistor (Th ) and the others (Th ), and 

I i 

normalizing with the variance of the first thermistor, 

c v = 1/T / Thi(t, x)Th-(t, x+ Axj) dt i=l 7 

1/T J (Th^t. x)) 2 dt 

where T is the time over which the cross variance was computed 
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(20 min). The spatial correlation was then plotted against separation 
distance between thermistors using thermistor one as a reference; 
curve fitting was done by linear interpolation. The unnormalized 
spatial correlation functions were then transformed to obtain wave 
number power spectra in the same manner as described for the 
temporal functions. 



/ 
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IV. RESULTS OF THE ANALYSIS 



A. TEMPORAL SCALES 

One of the most common ways of presenting scales of inhomo- 
geneities is to scale them by correlation times. In the present case, 
the correlation time is defined as the time required for the normalized 
autocorrelation function to fall to a value of 1/e (=0. 368). In water 
with a constant temperature gradient (such as the mixed layer), it 
might be expected that the shortest correlation times would be at the 
surface where turbulence and short period waves dominate the structure. 
With increasing depth where turbulent mixing is not so intense and 
small scale temperature gradients have begun to smooth out, the 
correlation times should lengthen out to reflect the dominance of the 
surface wave motion. In the vicinity of a thermocline having internal 
waves, the times should again increase to reflect the dominance of the 
longer period internal waves, unless the internal waves are creating 
turbulence, in which case the correlation times can shorten again. 

In the shallow water of this experiment, turbulent mixing (and there- 
fore shortening of the correlation times) may also have been caused 
by tidal currents. Consequently it is not surprising that the temporal 
correlations shown in table III show no dependence with depth. 

In addition, the vertical dimensions of the thermistor array were 
small with respect to the extent of the water masses involved. The 
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TABLE HI 



CORRELATION TIMES k SIGNAL VARIANCES 



Run 1 



Thermistor 


Time 

(sec) 


Var 

(°C 2 ) 

xlO' 5 


1 


11. 11 


301 


2 


11.72 


287 


3 


10.00 ' 


324 


4 


9.74 


337 


5 


8.81 


372 


6 


9.34 


471 


7 


4.37 


545 



Run 3 



Thermistor 


Time 

(sec) 


Var 

(°C 2 ) 

xlO" 5 


1 


3.28 


950 


2 


4. 10 


1021 


3 


3.69 


1230 


4 


3.89 


1611 


5 


10.73 


1722 


6 


5.40 


1295 


7 


5.09 


209 



Run 2 



Thermistor 


Time 


Var 




(sec) 


(°C 2 ) 

xlO" 5 


1 


9.74 


150 


2 


10. 80 


154 


3 


8. 59 


194 


4 


10.47 


178 


5 


8. 81 


153 


6 


9.75 


144 


7 


5.60 


91 


Run 4 




The rmistor 


Time 


Var 




(sec) 


(°C 2 ) 

xlO -5 


1 


14.40 


115 


2 


14.62 


112 


3 


12.88 


52 


4 


15.70 


115 


5 


17.26 


74 


6 


19. 54 


132 


7 


20.09 


65 



Thermistor 


Run 6 
Time 
(sec) 


Var, 

(°c 2 


1 


3.35 


xlO - 

.5660 


2 


3. 33 


5674 


3 


3. 16 


5449 


4 


3. 30 


3129 


5 


3.95 


4482 


6 


3.35 


11519 


7 


4. 60 


4784 
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array therefore never really spanned the different temperature regimes 
adequately to draw any conclusions. However, as will be seen, the 
possible influence of turbulent mixing at the therrmocline was notice- 
able on the traces of several thermistors on the edge of the array. 

Runs one and two were made with the vertical thermistor array 
at the top of the thermocline (Fig. 12, 13). The bottom thermistor 
of each run shows a shorter correlation time than those above. This 
reflects the possibility that turbulence was being generated at the ther- 
mocline in the vicinity of the lower thermistors (Fig. 1). It is interest- 
ing to note in passing the apparently strong surface wave influence on 
the thermocline in run one as represented in the signal from the bottom 
thermistor, number seven. The upper thermistor shows almost no 
surface wave influence. 

Run three, made in the thermocline, shows the most consistently 
short correlation times (Fig. 14) and tends to support the hypothesis 
that turbulence was being generated there. Discounting the correlation 
time for thermistor five which seems to be anomalously long, the two 
lower thermistors show longer correlation times than the others, 
indicating that perhaps the thermistor array extended below the 
thermocline outside of the active region. Run four, was made away 
from and above the thermocline and in nearly isothermal water 
(Fig. 15). It is the only one that shows long correlation times with 
respect to the other runs and shows a consistently increasing time with 
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0.368 



! 




Figure 12. Autocorrelation Functions Thermistors One and Seven Run One 
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0.368 




Figure 13. Autocorrelation Functions Thermistors One and Seven Run Two 



0.368 




Figure 14. Autocorrelation Functions Thermistors One and Seven Run Three 
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0.368 




Figure 15. Autocorrelation Functions Thermistors One and Seven Run Four 
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depth. This large difference between the correlation times for runs 
four and the others implies that it is the distance from the thermocline 
that causes this difference. Supporting evidence for this interpretation 
comes from the fact that coherence between the wave gauge and the 
top thermistor of runs one and four is different (Fig. 16). For run 
one (and two and three) the coherence at the principal wave frequency 
never exceeds 0.24 implying that the surface waves had little influence 
on microstructure movement. For run four, the coherence increases 
to 0.48 at the wave frequency implying that the surface waves had more 
effect on the microstructure movement in nearly isothermal conditions 
removed from the thermocline. The wave spectra for runs one and 
four are very similar in power and shape (Fig. 17); the increase in 
coherence is therefore not due to an increase in wave activity. The 
peaks in the coherence are narrow and the spectral energy shown by 
the thermistors of this run is greater at frequencies lower than the 
surface wave frequency implying that there was still energy being 
supplied from sources other than surface waves. These results in 
general support the findings of Whittemore [1973]. 

If the edges of the thermocline are assumed to be areas of 
increased turbulence, then thermistors in that area should show not 
only a reduced correlation time but an increase in variance over those 
far from the thermocline. The results of the experiment are not 
conclusive in this respect (Table III). The variance of the bottom 
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Frequency (Hz) 



Figure 16 




Frequency (Hz) 



Figure 17 



Comparison of Wave Spectra and Coherence 
Between Wave and Thermistor One Spectra 
Runs One and Four 
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thermistor for run one is nearly double that of the top thermistor 

(545 x 10 -5 °C vs 301 x 10 - ^ °C^) but for the second run it is smaller 

2 2 

by nearly half (91 x 10 ^ °C vs 150 x lO - ^ °C ). In run three, var- 

“5 2 

iances for the top six thermistors fall between 950 x 10 °C and 
5 2 

1722 x 10 °C but the bottom thermistor has a. variance of only 
-5 0 2 

209 x 10 C . An order of magnitude range seems to be a very 
large range, but a visual inspection of the strip chart shows that the 
bottom thermistor may have been outside the influence of the thermo- 
cline and the area of greatest turbulent activity. All the variances for 
run three which was made within the thermocline, are at least a factor 
of three greater than those of the other runs indicating that mixing 
may have been more intense than in the other runs. Variances for 
run four do not follow any consistent pattern and vary between 

2 *5 2 

65 x 10 ^ °C and 132 x 10 °C . No satisfactory explanation for this 
apparently random fluctuation could be found. 

Auto-spectra for selected thermistors were plotted in order to 
determine if any power law relationships existed between turbulence 
and frequency; signals with such a relationship would show up as 
straight lines on log-log plots. The plots (Fig. 18) show a tendency 
toward constant slope in the frequency range 0.2- 1.5 Hz, especially 
in the plot from run three thermistor one. The slopes are suggestive 
of turbulence in that the energy decreases with increasing frequency. 
However, the slopes are not consistent enough to support any generally 
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0 Thermistor 1 Run One 
X Thermistor 1 Run Two 
A Thermistor 1 Run Three 
© Thermistor 1 Run Four 
+ Thermistor 1 Run Six 



Figure 18. Autospectra for 
Temperature Signals, 

All Runs 
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Frequency (Hz) 

0.06 0.1 



accepted power law theory such as Kolmogorov^ inertial subrange law. 
The scales at which the data were taken were probably not sufficiently 
far removed from those of the internal waves to show purely turbulent 
characteristics . 

B. SPATIAL SCALES 

Thermistor signals were compared three ways in the spatial 
domain: spatial correlation, wave number spectra, and coherence. 

Spatial correlations for runs one and two are similar, with 
correlation lengths of 102 cm, reflecting their same general relation- 
ship with the thermocline. Run three shows a shorter correlation 
length of 33 cm; run four shows the longest correlation length of 
142 cm. Plots are shown in figure 19. 

The runs (one and two) made at the edge of the thermocline have 
longer correlation lengths than that of run three which was entirely 
within the thermocline, but shorter than that of run four which was 
above the thermocline. 

The wave number spectra are inconclusive. The spectrum for 
run one is figure 20. The energy from the very low wave number 
internal waves is apparently so great, even after filtering out the 
100 second period waves, that energy leaks into the area of higher 
wave number surface waves. 

For runs one through four and six, and for particular frequencies, 
coherency squared (or, simply, coherence) between thermistors was 
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Figure 19- Spatial Correlations for 
Temperature Signals, All Runs 
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Spectral Energy Density (xlO 
— , 1 — ■ 



o 

o 




Figure 20. Wave Number Spectrum Run One 
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plotted against thermistor spacing using thermistor one as a 
reference. Figure 21 is for 0. 008 Hz which may be considered as 
representative of the frequencies at which internal wave activity was 
found. Figure 22 is for 0.08 Hz where in general, maximum surface 
wave energy was found. Figure 23 is for 0. 29 Hz, an arbitrarily 
picked high frequency. At the lowest frequency, coherence for runs 
one, two, and four falls to 1/e (0.368) in approximately 90 cm. The 
coherence in run three falls to 0.368 in 17 cm. It can be seen again 
that the random temperature fluctuations caused by the turbulence 
associated with the internal wave have a drastic effect in reducing the 
temperature field coherence as they did in reducing the temporal and 
spatial correlations. 

At the surface wave frequency the picture is not so distinct. The 
run showing coherence of at least 0. 368 for the longest distance is, 
as expected from previous discussion, run four at 38 cm. Runs one 
and two show 0.368 or better coherence to 15 cm. Run three, in the 
thermocline and subject to an internal wave, shows coherence out to 
25 cm. 

At the high frequency chosen, the spacing at which coherence 
remains greater than 1/e is expectedly short. Any turbulence 
present will dominate this frequency. In general the runs made in 
the presence of internal waves and the turbulence they cause show loss 
of coherence in shorter distances than those less affected by internal 
waves. Loss of coherence occurs between 5 and 10 cm. 
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as a Function of Thermistor Spacing: 0.29 Hz 



C. COMPARISON OF THE VERTICAL AND HORIZONTAL SCALES 

The relation between the vertical and horizontal scale sizes can 
best be seen through a comparison of run three vertical array and run 
six horizontal array, since both of these runs display a great deal of 
internal wave activity. It should be kept in mind for this comparison 
that the thermocline was more intense for run three than for run six 
(0. 79°C/m vs 0. 44°C/m); also, the upper part of the vertical array in 
run three was subject to internal waves whereas a visual inspection of 
the trace for run six reveals that all thermistors were affected by 
internal waves for that run. Nonetheless, since internal waves can be 
major contributors to turbulence, the comparison will be made. 

A comparison of correlation times in table III shows that both runs 
have correlation times on the order of 3-4 sec. Variances for run 
six are, in general, consistent with each other, but they are about four 
times greater than those in run three. This may imply that mixing 
was more intense in run six. 

Spatial correlations for the two runs are remarkably similar, 
both showing correlation lengths of about 38 cm (Fig. 24). The 
coherence for run six at the internal wave frequency (0. 008 Hz) 
remains greater than 1/e over a much greater distance (58 cm) than 
for run three (17 cm) (Fig. 21). The implication is that internal 
waves are coherent over a longer distance horizontally than vertically. 
The coherence of the thermistors in run six at the surface wave 
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Figure 24. Autocorrelation Functions 
Thermistors One and Seven Run Six 
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frequency falls off surprisingly fast with thermistor spacing. A plot 
of the coherence for the wave gauge and thermistor five (almost 
directly under the wave gauge) shows one of the highest values obtained 
at the principal wave energy frequency during the experiment 
(0.40 at 0.08 Hz) (Fig. 25). This implies that coherence between the 
thermistors should be high at that frequency. This does not occur, 
possibly because the thermistor spacing brackets the lengths of the 
surface wave particle orbits at which maximum coherence occurs. 

All that shows is a very small peak (0. 26) at 96 cm. The coherence 
at 0.29 Hz is negligible and falls to 1/e at 5 cm. The fact that this 
value matches the value obtained from the runs with vertical arrays 
indicates the possibility of well developed turbulence at that frequency. 

D. INTERNAL WAVES 

Three aspects of internal waves could be seen from the data. 

First, when the thermistor array was vertical (run three), a wave 
passed which showed an instantaneous vertical picture of turbulence 
at the trough. Second, when the array was horizontal (run six), a 
trough passed giving a serial view of the wave energy being progres- 
sively degraded by turbulence. Third, also with the array horizontal 
(run six), the midpoint of an apparently large wave passed from 
which could be estimated the maximum probable dimensions of the 
wave and its phase speed. 



61 



Coherence 




Figure 25. Coherence of Signals from 
Wave Gauge and Thermistor Five Run Six 
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In the first case (Fig. 9) the temperature excursion across the 

wave is 2. 7 °C corresponding to a wave height of 3.4 m assuming a 

o 

constant temperature gradient of 0.79 C/m. Since the wave appears 
distinctly only in the traces from the top three thermistors of the array 
and their total spacing is only 15.2 cm, there are not sufficient phase 
shifts to compute phase speed. The trace from thermistor five shows 
the internal wave breaking down into turbulence; four distinct eddies of 
amplitude 0. 84 °C can be readily seen. The trace from thermistor 
seven shows the internal wave only just visible as a smooth diffuse 
undulation. 

The second view of an internal wave creating turbulence is in 
figure 10. The trough of the wave apparently shows as it passes the 
horizontal array. In the trace from thermistor one, the warm sector 
has already been subject to some mixing with colder water and shows 
considerable structure. The trace from the seventh thermistor shows 
the warm sector nearly completely destroyed as an entity; instead a 
group of four well formed eddies has taken its place. Phase relation- 
ships show this change occurred in about 12 sec. 

The third internal wave shown follows the one just discussed on 
the same run (Fig. 11). Since the array was horizontal, the thermistor 
traces again give a serial view of the wave passing. The traces are 

all similar in this case. The total temperature excursion of thermistor 

o 

one as the wave passes is about 2. 3 C; this corresponds to a wave 
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height of 5 m assuming a temperature gradient of 0. 44 °C/m in the 
thermocline. Phase relationships show a phase speed of 9 cm/sec. 

The warm sector is remarkably homogeneous and isothermal; tem- 
perature fluctuations are on the order of 0. 14 °C whereas the fluc- 
tuations in the surrounding cold sections are on the order of 0.42°C. 
This homogeneous warm water apparently comes from the surface 
mixed layer whereas the cold water comes from the vicinity of the 
thermocline where turbulence and microstructure are to be expected. 
The sudden onset and cutoff of intense temperature fluctuations (on 
the order of 1.4°C in three seconds) on the boundaries of the wave 
are indicative of intense turbulent mixing. However, the mixing does 
not appear to be sufficiently advanced to have caused any degradation 
of the warm sector. This wave possibly is not subject to as intense 
a density gradient as the preceding wave due to the turbulent mixing 
caused by that preceding wave. BT casts were not made with sufficient 
frequency to show such short term fluctuations to determine this. 
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V. CONCLUSIONS 



1. The temporal scales measured in the horizontal were greatly- 
affected by turbulence. In the presence of turbulence, scales were of 
the order 4 seconds with a standard deviation of 0. 8 seconds. Outside 
the influence of the internal waves scales were on the order of 10 
seconds with a standard deviation of 4 seconds. In the presence of 
turbulence, correlation times were similar for both vertical'and 
horizontal scales. 

2. The spatial scales outside the influence of turbulence were of 
the order of 118 cm with a standard deviation of 21. 5 cm. Under the 
influence of turbulence, they shortened to 33 cm. Vertical and hor- 
izontal scales in turbulent water were again similar. 

3. The principal cause of turbulence appears to have been 
internal waves. 

4. The data show the complexity of the temperature field off 

the NUC tower; they give some idea of the interplay between turbulence 
and internal waves during conditions of light winds, small waves and 
swell. 
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VI. RECOMMENDATIONS FOR FURTHER RESEARCH 



1. It became apparent during this analysis that much firmer 
conclusions could have been drawn if more frequent BTs had been 
taken or a more continuous monitoring of the vertical temperature 
structure could have been made. A vertical thermistor array nor- 
mally operating at the tower was not available during this experiment. 
Its use would have greatly aided interpretation of the data. Specif- 
ically, in some of the runs, events occurred at the top thermistors 
before the lower ones. This implies an internal wave cresting and 
beginning to break. Without the vertical profiles, this idea remains 
little more than speculation. 

2. Current measurements at several depths would have given 

a measure of velocity shear and aided interpretations of the turbulent 
and internal wave fluctuations. 
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